Testing isotropy of the universe using the Ramsey resonance technique on ultracold 

neutron spins 
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Physics at the Planck scale could be revealed by looking for tiny violations of fundamental sym- 
metries in low energy experiments. In 2008, a sensitive test of the isotropy of the Universe using has 
been performed with stored ultracold neutrons (UCN) , this is the first clock-comparison experiment 
performed with free neutrons. During several days we monitored the Larmor frequency of neutron 
spins in a weak magnetic field using the Ramsey resonance technique. An non-zero cosmic axial 
field, violating rotational symmetry, would induce a daily variation of the precession frequency. Our 
null result constitutes one of the most stringent tests of Lorentz invariance to date. 



I. INTRODUCTION 

The concept of symmetry plays a central role in our 
current understanding of physics at its most fundamen- 
tal level. The invariance of physical laws under rotation 
symmetry is perhaps the most basic example, and for 
such it is worthwhile to test its validity with high accu- 
racy. Moreover, it has been argued that a breakdown of 
rotation symmetry, or more generally Lorentz symmetry, 
could arise from quantum gravity effects. The high level 
of precision reached by modern tests of rotation symme- 
try opens the possibility to observe tiny manifestations 
of the physics at the Planck scale. Here we describe an 
experiment performed recently with ultracold neutrons 
based on spin dynamics Picture a neutron sitting in 
an environment free of magnetic field (and electric field) . 
One consequence of rotation invariance is the fact that 
the neutron spin should be stationary, since there is no 
direction around which it could turn around. If the neu- 
tron spin were rotating around a direction b at a nonzero 
frequency -\b (H is the Planck constant), we would ob- 
serve a breakdown of the rotation invariance. Such an 
effect is parametrized by the exotic potential V = o~ ■ b, 
where a are the Pauli matrices acting on the neutron spin 
states. The vector b, so-called cosmic spin anisotropy 
field, defines a priviliged direction in the universe. Now 
if a neutron is subjected to both a static magnetic field B 



and the new field b, its spin will precess at the modified 
Larmor frequency f n , which to first order in b is given 
by: 

f n= 7n. B + ±. h .*. (1) 

If we keep the magnetic field B static in the laboratory 
frame fixed on Earth, say vertical, then the scalar prod- 
uct b • B will be daily modulated as the Earth is rotating. 
Searching for a daily modulation in the Larmor frequency 
in a static magnetic field probes the component of b 
orthogonal to the Earth's rotation axis. 

The presented result is in fact a clock comparison ex- 
periment in that the neutron Larmor frequency is com- 
pared to the precession frequency of a co-habiting 199 Hg 
magnetometer. In the last section of this article we will 
discuss our result in regard with the numerous previ- 
ous clock comparison experiments, in particular those 
performed with 199 Hg. The first clock comparison ex- 
periments were performed by Hughes Q and Drever 
They were motivated by testing Mach's principle, aiming 
at detecting a dependence of the Zccman energy of nuclei 
as a function of the relative orientation of the direction 
of the magnetic field and the direction to the galactic 
center. Nowadays these searches follow the development 
of a general formalism for Lorentz violation called the 
Standard Model Extension (SME) Q that parametrizes 
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all Lorentz-violating effects possibly induced by the yet 
unknown Planck scale physics. The SME contains a num- 
ber of terms that violate local Lorentz invariance, among 
which is the cosmic axial field b. 

Experimental setups to resolve changes in neutron Lar- 
mor precession frequency have been in constant develop- 
ment over the past 50 years with the goal of detecting a 
non zero neutron electric dipole moment. For this test of 
rotation symmetry we used the most sentive one of such 
apparatus available today, namely the OILL apparatus. 
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FIG. 1: Sketch of the OILL apparatus. 
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FIG. 2: Number of UCN counts in the detector as a function 
of the storage time in the OILL chamber. 



can be stored in material bottles for a very long time. 
In a typical cycle the UCNs are guided from the source 
to a storage volume (a cylinder, 12 cm height and 47 cm 
diameter) , filling the volume for about 40 s before closing 
the UCN valve. On the way to the storage volume they 
are polarized when passing through a magnetized ferro- 
magnetic foil. After the completion of the Ramsey pro- 
cedure described below, the UCN valve is opened again 
and the neutrons fall down into a He detector where they 
are counted. On the way to the detector they pass again 
through the magnetized foil that serves as a spin anal- 
yser. Figure [2] shows the number of UCN counts as a 
function of the duration of the storage time. From the 
nearly exponential decrease of the neutron count rate we 
extract the UCN storage time of Tq 200 s in the OILL 
chamber. The storage curve departs from a simple ex- 
ponential behaviour because different parts of the UCN 
energy spectrum have different storage times. 



II. THE OILL SPECTROMETER 

The OILL spectrometer has been built and operated 
at the PF2 beamline at ILL by the Sussex/RAL/ILL col- 
laboration and it holds the world record limit on the neu- 
tron electric dipole moment [5j . It is now operated by an 
European collaboration in view of a more sensitive mea- 
surement Q at the Paul Scherrer Institut where it will 
profit from a more intense ultracold neutron source. The 
data reported in this article has been collected at the ILL 
in 2008, the year preceding the move of the apparatus. 
Figure [T] shows schematically the OILL spectrometer. 



A. Storing polarized ultracold neutrons 

The OILL apparatus uses ultracold neutrons, neutrons 
with kinetic energy of about 100 neV. Those neutrons are 
reflected at all angles of incidence by most materials and 



B. Measuring the neutron Larmor frequency 

The chamber is exposed to a static vertical magnetic 
field of Bq = 1 /iT, corresponding to a Larmor frequency 
of fn ~ 30 Hz. During the storage of polarized UCN the 
Ramsey method of separated oscillatory fields is applied, 
in order to measure /„ accurately. An initial oscillating 
horizontal field pulse of frequency /rf is applied for 2 s. 
It flips the neutron spin by 7r/2. Then the UCN spins pre- 
cess freely around the Bq field, for typically T = 100 s. 
A second 7r/2 pulse, in phase with the first pulse, is then 
applied. The Ramsey procedure is resonant. It flips the 
neutron spin by tt only when /rf = fn- Figure [3] shows 
a measured Ramsey resonance: the vertical component 
of UCN spins as a function of the pulse frequency /rf- 
To optimize the sensitivity to a change in the Larmor fre- 
quency we operate at /rf corresponding to equal spin up 
and spin down counts in the sides of the central Ramsey 
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FIG. 5: Variations of the R ratio around its average. For 
clarity the data are folded modulo 23.9345 hours and binned 
every half hour. 



FIG. 3: Ramsey resonance measured with OILL. Red: counts 
of the spin up neutrons, black: counts of the spin down neu- 
trons. 
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FIG. 4: Each point is a measurement cycle of the neutron 
Larmor frequency, uncorrected (dots) and corrected for the B 
field fluctuations with the comagnetometer (triangles). 



fringe where the slope is the steepest. One can show that 
each cycle gives a measurement of the Larmor frequency 
f n with a statistical accuracy of: 
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where N is the number of neutron counts, and a is the 
polarization of the neutrons at the end of the storage 
period. In operating conditions this was of the order of 
(Tf n rj 30 /iHz. 



C. Controlling the magnetic field 

To resolve a daily modulation in the neutron Lar- 
mor frequency, one has to make sure that the magnetic 
field Bq is not fluctuating too much. The control of the 
magnetic field is the most critical issue for nEDM spec- 
trometers and involves a great deal of shielding, shaping 



and monitoring of the magnetic field. First, the storage 
chamber is surrounded by a four layer mumetal mag- 
netic shield that protects the experiment against varia- 
tions of the external field. The shielding factor against 
slow changes of the outside vertical field was measured 
to be about 10000. But this is far from being enough 
to lower the inner field fluctuations below the statistical 
accuracy ([2]) of the neutron cycles, as is illustrated in fig. 
2] To gain control over the residual magnetic field fluctu- 
ations, a unique feature of the OILL spectrometer is the 
mercury comagnetometer [8|. Within the neutron stor- 
age chamber, nuclear spin-polarized 199 Hg atoms precess 
in the same magnetic field as the neutrons. The Larmor 
frequency /n g ~ 8 Hz of the mercury atoms is measured 
optically, by recording the oscillation of the transmission 
of polarized resonant light. The mercury comagnetome- 
ter provides a field average measurement for each cycle, 
with the same time average as for the neutrons, at an 
accuracy of tr/ Hg « 1 /iHz. The performance of such 
a method is illustrated in fig. [4] where the neutron fre- 
quency is corrected for changes of the magnetic field using 
the comagnetometer. The fluctuations of the corrected 
neutron frequency are in agreement with the purely sta- 
tistical precision ((2]). 



III. RESULTS AND STATISTICAL ANALYSIS 

The fact that the fluctuations of the magnetic field 
are suppressed below the statistical sensitivity was used 
to search for a daily modulation of the neutron Larmor 
frequency. For each measurement cycle we form the ratio 
of the neutron to mercury Larmor frequency R = fn/fng- 
This ratio is free from magnetic field fluctuations. Then 
we parametrize the signal for the daily modulation by 
the amplitude A and the phase 4>: 
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where ft = 2ir/23. 9345 his the sidereal angular frequency, 
7m 7Hg are the gyromagnetic ratio of the neutron and the 
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FIG. 6: Bayesian probability density for the amplitude A. 



mercury atoms, and SR is a constant systematic shift that 
can occur for example due to magnetic field gradients, 
discussed in the next section. 

In 2008 two series of runs have been conducted. The 
first series, recorded in April-May 2008 contains in to- 
tal 4.6 days of data in three different continuous runs. 
This series was analysed [l[ and no daily modulation 
was found in the ratio R. An upper limit on the max- 
imum amplitude of the modulation was reported A < 
0.58 x 10 -6 (95 % C.L.) using a frequentist statistical 
method. In magnetic field units this would correspond 
to Ba — ^4Bo7Hg/7n = 150 fT. A second series was 
recorded in December 2008, consisting of three contin- 
uous runs of 5.6 days in total. In this series an electric 
field was applied, with reversed direction every 2 hours 
for a nEDM measurement. Since the frequency of the E 
field reversal is large compared to the sidereal frequency 
the electric field does not affect this analysis even if the 
neutron EDM were enormous. 

In this article we present a combined analysis of both 
series, searching for a daily modulation. For each run 
the average value of R has been subtracted, resulting in a 
subtracted value AR for each cycle, with a standard error 
or. Fig. [5] shows an overview of all data, folded modulo 
23.9345 h. The error bars indicate combined statistical 
errors of the neutron and the Hg frequency, dominated 
by the former one. The whole dataset is compatible with 
a signal of null amplitude (reduced x„ u n = 0.90). 

A Bayesian analysis was applied to the data to search 
for a time variation R(t), Eq. ([3]). First, the following 
Chi squared function is established: 

*WHE ( ^-^ + *) ) 3 , (4) 




0.4m 



0.5m 



FIG. 7: Vertical cut through the cylindrical storage chamber 
for UCN and 199 Hg. Schematically indicated are the « 70 mm 
diameter Cs vapor filled bulbs and their mounts. The scalar 
Cs magnetometer measures the magnitude of B found at the 
center of the spherical bulb. 

function: 

L(A) = 1 J 2 J exp(- x 2 (A, </,)) d<j>, (5) 

where N is a normalisation coefficient. This function is 
plotted in Fig. [6] from which we extract the following 
bound: 

A < 0.35 x I0~ 6 95 % C.L. (stat) (6) 

It improves the bound obtained with the first series alone 
as expected. 

IV. POSSIBLE SYSTEMATIC EFFECT 

A systematic effect could arise due to imperfect mon- 
itoring of the magnetic field by the comagnetometer. If 
the source of this imperfection is daily modulated it could 
in principle hide a true cosmic signal. The only known 
relevant effect is the so-called gravitational shift. The 
mercury atoms form a gas at room temperature and fill 
the storage volume with a uniform density. On the con- 
trary, the UCN gas is affected by gravity and the UCN 
density is significantly higher at the bottom of the stor- 
age chamber than at the top. This results in a difference 
in center of mass height h of the two species, that in 
turn results in a shift of the ratio R in the presence of a 
vertical gradient dB/dz of the magnetic field: 

7H g V B dz J 



where the sum runs over the data cycles. The poste- 
rior probability density for A is given by the likelihood 



A first estimate of h can be obtained from the kinetic 
theory of gases. The maximum energy of stored UCN in 
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FIG. 8: Correlation between the vertical field gradient mea- 
sured by the Cesium magnetometers and the R = f n /fn e FIG. 9: Fluctuations of the magnetic field gradient measured 
ratio, for magnetic field pointing downwards (blue) and up- by the Cesium magnetometers, 
wards (red). 



the chamber is the Fermi potential of the weakest wall, 
i.e. Vf = 95 ncV for the quartz insulator ring. We 
can assign a temperature Ttjcn to the UCN gas by the 
matching of the mean energy: I^bTtjcn = § Vf, that 
is, Ttjcn — 0.5 mK. Assuming "thermal" equilibrium for 
the UCN gas the distribution of height can be derived 
together with the center of mass offset: : 



f{z)dz 
h 



= exp 



J H (H/2-z)f(z)dz 
f Q f{z)dz 



mgz 



J dz, 



= 2.9 mm, 



(8) 
(9) 



where H = 12 cm is the height of the storage cham- 
ber. This is only an estimate since the UCN energies 
are certainly not distributed according to the Maxwell- 
Boltzmann distribution. 

A dedicated experiment was performed to directly ob- 
serve the gravitational shift (0 using four Cesium mag- 
netometers Q running in parallel with the neutron and 
mercury systems. These magnetometers were placed as 
described in fig. [7J two above the storage chamber and 
two below, allowing to measure the magnetic field gra- 
dient dB/dz. The magnetic field gradient could be ad- 
justed using correction coils at the top and at the bot- 
tom of the vacuum chamber. Figure [8] shows the correla- 
tion between the gradient extracted from the Cs magne- 
tometers and the ratio R. Each point corresponds to a 
different magnetic field configuration, i.e. different cur- 
rents in the correction coils. The correlation is that ex- 
pected from |7|). with a measured center of gravity offset 
of h = 2.3 ±0.1 mm. This is in fairly good agreement 
with the naive estimate ©. 

A daily modulation of the magnetic field gradient 
would induce a false modulation of the R ratio through 
eq. ([7]). During part of the datataking, the gradients 
were measured online with the Cesium magnetometers, 



as shown in fig. [9j From this measurement we extracted 
the daily modulated Fourier component of this signal, 
with an amplitude of 30 pT/m. The statistical analysis 
presented in the previous section was repeated, this time 
with a modulation of amplitude 0.3 ppm superimposed 
to the data. The worse case scenario (the worse phase 
of the superimposed oscillation) leads to the corrected 
bound: 



A < 0.4 x 10~ 6 95 % C.L. (stat + syst). 



(10) 



The systematic effect just discussed is still negligible, but 
would become an issue for future experiments at the Paul 
Scherrer Institut with improved statistics. 



V. DISCUSSION 

The result (jTUJ) can be interpreted in terms of a limit on 
the cosmic spin anisotropy field b for the free neutron. 
The amplitude A is related to the component b± of b 
orthogonal to the Earth rotation axis: 



A = b 



cos(A) 



(11) 



where A is the latitude of the experiment. Thus we ex- 
tract a limit on the small energy scale associated to the 
breakdown of rotation symmetry: 



6 ± <lxl0~ 20 eV 95% C.L. 



(12) 



This limit constitutes a improvement by a factor of two 
compared to our previous analysis [l[ based on only part 
of the data (b± < 2 x 10~ 20 eV). The new dataset added 
(December 2008) to the analysis was taken 7 months 
apart from the first dataset (April-May 2008). This is 
ideal to cancel any day-night fluctuations because the 
day-night phase and the sidereal phase are in n phase 
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opposition after 6 months. In this respect, the new com- 
bined analysis (|12p is free of unknown day-night influ- 
ence. 

Table U compares the present result (fl"2"j) to existing 
limits on other particles. The result reported here is the 
best limit for the free neutron. It is complementary to 
the more precise atomic experiments [ToL 111! that can be 
interpreted as limits concerning bound neutrons inside 
nuclei. Our result using UCN is free from model depen- 
dent nuclear corrections and related possible suppression 
effects. In our analysis we have assumed that the Larmor 
frequency of the mercury atoms was not daily fluctuating 
due to the cosmic spin anisotrop y fi eld. This assumption 
is justified by the existing limit [l(| for 199 Hg. 



TABLE I: Results of more restricting upper limits (at 
95% C.L.) on b±(e), b±(N), b±(p), b±(n), b±(n) the couplings 
between a cosmic spin anisotropy field and different particles. 



Reference 



System 



Particle 



bx [eV] 



Berglund et al., [10J Hg & Cs bound neutron 9 x 10 



electron 



2 x 10 



-20 



Bear et al., fill 
Phillips et a/., [12] 
Heckel et al., [131 
Bennet et al., [14] 

This analysis 



Xe & He bound neutron 2 x 10~ 22 
H proton 4 x 10" 18 

e electron 7 x 10~ 22 

(j, positive muon 2 x 10~ 15 
negative muon 3 x 10~ 

n & Hg free neutron 1 x 10~ 20 



We have shown that if a neutron would sit in a zero 
magnetic field environment, the period for the rotation 
of the spin around a priviledged direction in the universe, 
if any, would certainly be longer than 2 days. 
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